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We study the renormalization group flow of Z 2 -invariant supersymmetric and non-supersymmetric 
scalar models in the local potential approximation using functional renormalization group methods. 

We focus our attention on the fixed points of the renormalization group flow of these models, which 
emerge as scaling solutions. In two dimensions these solutions are interpreted as the minimal (su¬ 
persymmetric) models of conformal field theory, while in three dimensions they are manifestations 
of the Wilson-Fisher universality class and its supersymmetric counterpart. We also study the ana¬ 
lytically continued flow in fractal dimensions between 2 and 4 and determine the critical dimensions 
for which irrelevant operators become relevant and change the universality class of the scaling so¬ 
lution. We include novel analytic and numerical investigations of the properties that determine the 
occurrence of the scaling solutions within the method. For each solution we offer new techniques to 
compute the spectrum of the deformations and obtain the corresponding critical exponents. 


I. INTRODUCTION 

The birth of functional renormalization group methods 
is tightly connected to the investigation of scaling prop¬ 
erties of the effective potential in statistical field theories 
with Z 2 or O (N) symmetries jU [2j, and it represents its 
first non-trivial application [3H5]. Over the years the in¬ 
terest in properties that can be extracted from the renor¬ 
malization group flow of the potential in a local approx¬ 
imation has been revitalized at various stages IS El- 

The effective potential of a Z 2 -symmetric scalar field 
theory has been studied with various techniques to un¬ 
veil the underlying universality classes, which include the 
Ising model. The critical exponents have been deter¬ 
mined to high accuracy within polynomial expansions W 
E] and non-polynomial solutions have been studied nu¬ 
merically by applying scaling solutions methods mm- 
Recently, even global solutions have been constructed 
with very high numerical precision |15| . 

In the past few years, the scaling solutions approach 
was adopted to analytically continue various parameters 
which enter the renormalization group flow, including the 
dimensionality of space |T6j, the number N of fields pn- 
[2P] or both [2Tl. The analytic continuation of these pa¬ 
rameters offers new insights into properties of the un¬ 
derlying statistical models, such as their emergence on 
spaces of fractional dimensions, the occurrence of the 
Mermin-Wagner theorem and constraints on the exis¬ 
tence of the various multicritical models. 

In this paper we extend the scaling solution approach 
in three prominent ways. First, we will carefully lay down 
and discuss the theoretical mechanism with which a scal¬ 
ing solution occurs for a specific statistical model, which 
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has often been overlooked in previous works. This occur¬ 
rence will be illustrated mostly through the example of 
the three dimensional Ising model. Second, we will in¬ 
troduce two methods to compute the critical exponents 
of a model: a shooting approach (which has been used 
previously in [221424] f to calculate critical exponents in 
three and two dimensions) and the novel pseudo-spectral 
method with a SLAC derivative f!o\f27\ . These two meth¬ 
ods treat the boundary conditions for the deformations 
of a critical potential in different ways; thus, they are 
complementary for the approach as will be shown by ap¬ 
plying them to the known Ising example. Finally, we 
systematically apply the scaling solution approach to the 
supersymmetric N = 1 Wess-Zumino model in both two 
El [301 and three dimensions |31H33j . for which we com¬ 
pute various critical exponents for its universality classes. 
We also study the supersymmetric model on fractional 
dimensions introducing a novel analytic continuation. 

The paper is organized as follows: The remaining part 
of Section [I] introduces the local potential approximation 
and the functional renormalization approach, Sect. |TT] in¬ 
troduces the scaling solution approach and applies it to 
various universality classes in two and three dimensions, 
Sect. DU introduces the Wess-Zumino model and shows 
the application of the scaling solution method to the su¬ 
persymmetric analogues of the universality classes of the 
previous Section, Sect. |IV| shows how the universality 
classes emerge at certain critical dimensions when vary¬ 
ing the dimensionality of the system, and finally Sect. [V] 
briefly discusses our results. The appendixes are dedi¬ 
cated to further clarifying some properties of the renor¬ 
malization group flow and the methods used in this work. 


A. Local potential approximation 

For clarity purposes we introduce all the basic concepts 
used in the present paper adopting the illustrative exam- 
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pie of a simple Z 2 -invariant scalar field on This intro¬ 
duction can and will be straightforwardly generalized to 
the inclusion of supersymmetry later in the correspond¬ 
ing Sections. Readers that are experts in the techniques 
of Wilsonian and functional renormalization are welcome 
to skip the remainder of the present section. 

The local potential approximation (LPA) is a non- 
perturbative truncation of the effective action of a system 
which includes all the relevant operators in a low energy 
derivative expansion 1 13 t l.'il l . Within this approximation, 
higher derivative operators are neglected but all local 
interactions are included in a non-perturbative fashion. 
The LPA of a simple Z 2 -invariant scalar p in d dimen¬ 
sions thus truncates the effective action of the system 
to 


m = J d d x\^d p pd> 1 p + V{p)^ , ( 1 ) 

where we included a wavefunction renormalization con¬ 
stant Z for the field p and the local potential V(p). More 
precisely, the truncation 0 represents an extension of 
the original LPA of a simple scalar that goes under the 
name of improved LPA (a.k.a. LPA’), as it includes a 
non-trivial wavefunction renormalization for p and thus 
allows the field to have a non-canonical scaling in the 
guise of an anomalous dimension for the renormalized 
field ipii = Z 1 / 2 p. 

B. Functional Renormalization 

We study the renormalization group flow of the scalar 
models using functional renormalization group (FRG) 
methods. In this approach the path-integral of the sys¬ 
tem is constructed as 

e -r*M = 

where S[x] is a bare action for the microscopic field \ 
and A,S’ fe [x] is a cutoff that suppresses the propagation 
of the infrared modes relative to a given reference scale 
k. The infrared cutoff should be quadratic in the field 

A5fc[x] = \ j d d xd d x' x(x)K k (x - x')x(x') , (3) 

and its kernel 7 Z k suppresses the propagation of the in¬ 
frared modes by directly modifying the bare propagator 
of the theory. The dependence of 0 on the kernel 7 Z k 
encodes the scheme-dependence of the method; thus, dif¬ 
ferent kernels implement different coarse-grainings. How¬ 
ever, the kernel has to satisfy in the limit k —► 0 the 
property TZ k ^o = 0 , which ensures that 0 reduces to 
the standard path-integral of the (full) effective action 
T[p] of the system when the reference scale k is taken to 
zero. 

The functional T/ c [p\ appearing on the rhs of ([2]) is 
known as the effective average action (EAA) of the 




( 2 ) 


model. By construction it reduces to the full effective 
action T^oM = r[<p] in the infrared limit k —> 0. Fur¬ 
thermore, for sufficiently regular cutoff kernels it can 
be shown that an ultraviolet limit can be constructed 
T^aM = <Sa [</>], in which A is an ultraviolet cutoff and 
SaM i s the (possibly divergent) bare action of the model 
at the UV-scale A. While the path integral 0 is gen¬ 
erally not free of ultraviolet divergences and, in princi¬ 
ple, requires renormalization, it is possible to show that 
the flow of F* [<p] with respect to the scale k is finite. 
Therefore, rather than achieving the renormalization of 
the theory through the inclusion of counterterms, it is 
customary to use (J2| to derive the exact and finite FRG 
equation 

kd k r k [<p\ = ^Tr(r ( * ) [<p\ + n k s ) 1 kd k n k , (4) 

and use it to investigate the renormalization of the model. 

The flow encoded by Q determines a velocity held in 
the space of all functionals of the Helds entering a model 
under investigation. The LPA can be applied easily by 
projecting both sides of Q to a fc-dependent generaliza¬ 
tion of ([!]) which reads 

r fe [<d = J d d z j^d /i <pd'V + R fc (v?)j • ( 5 ) 


The anomalous dimension of the Held is 77 = —kd k Z k /Z k . 
In order to read off the renormalization group flow of the 
couplings in ([5]) we choose the prescriptions 


kd k V k (p) = ^fcd fe r fe [p]|^ =const , 


J_ fc& A 

Z k dp 2 5p p 5p_ p 




( 6 ) 


which must be applied to the right hand side of Q . Var¬ 
ious remarks are in order to explain the meaning of 0: 
The flow of the local potential is defined as the (normal¬ 
ized) flow of the EAA in the limit of constant field con¬ 
tent, as is usual in quantum field theoretical applications. 
The field p p represents the Fourier transform of p(x) with 
respect to the momentum p; thus, the anomalous dimen¬ 
sion 77 is extracted quite naturally from the coefficient of 
the p 2 term in the two-point function belonging to T k [p]. 
Since the truncation ([5]) includes only a field-independent 
wavefunction renormalization constant, it is necessary for 
consistency to project the two-point function of r*,[<p] to 
a distinguished field configuration po- In applications po 
is chosen in a physically meaningful way as the minimum 
of the effective potential defined by V£(po) = 0. The con¬ 
figuration po thus plays the role of the order parameter 
of the Z 2 -symmetry po = (p), and the possibility that 
it has non-zero values represents a clear separation from 
standard perturbation theory. 

Let us now define the modified momentum-space prop¬ 
agator at constant field 


g k (q 2 ) = ( Z k q 2 + V£(p) + 7 Z k (q 2 )) 1 , (7) 
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where lZ k (q 2 * ) is (27r) d times the Fourier transform of 
the cutoff kernel in coordinate space. With the modi¬ 
fied propagator we can give a compact evaluation of i> 
based on a general cutoff 

kd k V k (ip ) = 1 d J d d qG k kd k U k 

” = -W0T J iiq + * 1 ®) „ 

°( 8 ) 

Here, the prime denotes the derivative with respect to 
the argument q 2 , whereas elsewhere it is the derivative 
with respect to the (constant) field. 


II. SCALAR FIELD THEORY 


We can give explicit expressions for the flow equation 
for the effective potential and the anomalous dimension 
after specifying the cutoff function. While, in principle, 
any choice compatible with the properties of a suitable 
cutoff is possible, for the purpose of this paper it is useful 
to choose one that allows for an explicit calculation of 
the flow. We choose the popular optimized cutoff and 
explicitly factor out the wavefunction renormalization of 
the field, so that in momentum space the kernel of ([3| 
reads 

n k (q 2 ) = Z k (k 2 - q 2 )e(k 2 - q 2 ) . (9) 

The optimized cutoff is a distribution and therefore the 
momentum space integrals <[ 8 j) must be treated with 
careQ 

Let us also introduce the dimensionless renormalized 
field tp R = Z k 1 ^ 2 /c( 2 _d )/ 2 y;, which is the natural argu¬ 
ment of the dimensionless renormalized potential 

MPr) = k~ d V k {ip). (10) 

Together and ffc(<p R ) are the appropriate variables to 
investigate the scaling properties of the system, and in 
fact we shall use them for the remainder of this section. 
For notational simplicity we also drop both the labels 
and the overline notation, so the explicit evaluation of 


1 The presence of the derivatives with respect to q 2 in § requires 
a normalization for the integral of the product of the Heaviside 
step function 0(y — yo ) and its first derivative 5(y — yo). Given 
a regular test function g(y) we normalize 

J dy g{y) 8{y - yo) 6{y - y 0 ) = 1/2 9(1/0). 

This normalization is consistent with the limit of the same inte¬ 
gral evaluated for any sequence of smooth functions approaching 

0(y-yo)- 


the flow in the new variables takes the form 


kd k v{ip) = S VtTI [v,v';<p] + T V:v [v "], 

Sy, v [v, v’\<p] = -dv(<p) + - — Y^-<pv'{<p ), 
1 _ 3 _ 

„//l _ „ d+ 2 


-r \ _ a-\-z 

1 - Cd i + u"M ’ 


(ii) 


in which we defined the dimension-dependent geometric 
constant c^ 1 = (47r) d ' /2 r(l + d/2). We split the flow of 
the potential in two parts with different origins. The 
functional v'\ ip] represents the scaling part of the 

flow and it originates from the change of dimensionful 
variables to dimensionless ones. Its form is independent 
from the cutoff choice § and therefore from the scheme 
of computation. The functional J- V , v [v”], instead, en¬ 
codes the nontrivial features of the flow and depends on 
the cutoff. 

We now compute from ([ 8 j the anomalous dimension 
explicitly to obtain 


77 Cd {1 + v" (ip 0 )) 4 


( 12 ) 


where again we assume ipo to be the expectation value of 
the order parameter of the system. For later convenience 
we also give the result of an alternative computation of 
the anomalous dimension, which comes as the limit N —> 
1 of the anomalous dimension of the Goldstone modes of 
the 0(/V)-model 


?/ = 2 c d 


v"(y 0 ) 2 

¥>o(l + v"{tp 0 )) 2 ' 


(13) 


This definition has often and successfully been used in 
place of 77 on the basis of physical arguments. Thus, we 
include it in our analysis and refer to it as Goldstone 
anomalous dimension. The reason why 77 and r) 1 are dif¬ 
ferent is not because the limit N —> 1 of the 0(A)-model 
is discontinuous, but rather because 77 represents the lim¬ 
iting value of the anomalous dimension of the massive 
Higgs field in the broken O(TV) phase. 


A. Scaling solutions 


Let us now discuss some general property of the flow of 
the potential The computation can be manifestly 

analytically continued to any non-integer dimensionality, 
and for the present work we are mostly concerned with 
the interval 2 < d < 4 which is bounded by the lower 
and upper critical dimensions of the ip 4 model. Note that 
the flow of the potential depends parametrically on the 
anomalous dimension given in (12) or (13). Nevertheless, 


for 0 < 77 < 1 we can extract some general properties of 
the flow that are insensitive to the values of 77 . The flow 
(11) is a PDE representing the evolution of the potential 


with respect to the scale k. The fc-independent (station¬ 
ary) solutions are the critical solutions (fixed points) of 
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the flow, and they correspond to critical models of the 
corresponding field theory. This means that critical so¬ 
lutions solve the non-linear ODE 

S v ,ri[v,v';tp] +Tv,n[v"] = 0, (14) 

which is called a fixed point equation. The function J~ v ,r) 
can be explicitly inverted to cast the ODE in a standard 
form 


v"(tp) = -F^iSv^v'wW . (15) 

While the explicit form of this ODE depends on the cut¬ 
off, it is easy to prove that for any cutoff 

Tv,r,[x\ ~ 1 /s, for x —> oo , (16) 


which implies 

v"(p) ~ -l/S VtV [v,v'-,ip ], when S v , v ~ 0. (17) 

In other words the second derivative of the potential is 
expected to diverge at points p at which the scaling part 
becomes zero. These points can be determined implicitly 
by the equation S v = 0 as 


p = 


2d v(ip) 
d — 2 + r) v'(p ) 


(18) 


It is possible to show that any attempt to integrate the 
fixed point equation (15) starting from initial conditions 


at p = 0 will stop (with probability one) at a singularity 
located at p at which v"{p) diverges, but v(p) and v'(p) 
stay finite [22j . In addition, if p > 0 then v(p) and v'(p) 
have the same sign; similarly, if p < 0 they have opposite 
signs. 

Ideally, we are interested in global analytic solutions 


of (15), which thus avoid the appearance of singularities 
such as (18) for finite values of the field p, and we hope 


that these global solutions can be obtained by numeri¬ 
cally integrating the ODE (15) with fine tuned boundary 


conditions. In practice, it turns out that in the process 
of integrating (15) singularities are the rule and not the 


exception. To investigate the singular solutions we first 
choose a suitable parametrization. Since (15) is a second 


order ODE, any solution can be characterized by the two 
initial conditions u(0) and i/(0). One of the conditions is, 
however, fixed by the requirement of Z 2 -invariance to be 
v'(0) = 0 , thus ensuring the even parity of the solutions. 
It proves convenient to parametrize the second condition 
using the value u"(0) rather than u(0). We thus choose, 
as an initial condition i/ 7 ( 0 ) = cr, which implies that 


1 


u(0) = —-II— „ 

Cdd \ d + 2 


V 


1 


1 + <J ’ 


(19) 


which makes explicit the fact that cr = — 1 is a singu¬ 
lar point for the ODE. This comes from the more gen¬ 
eral property that, for any cutoff, the flow is well de¬ 
fined whenever v"(p) > —1. With the exception of the 


trivial case cr = 0, it is possible to observe numerically 
that almost all solutions that can be considered within 
2 < d < 4 diverge at a singular point p defined by (18) 
whose position depends on a. We therefore denote p„ 
the location of the singularity as a function of cr. 

The second derivative v"(p) diverges at field values for 
which the scaling part S Vir) vanishes. At any maximum 
or minimum of the potential with positive v the scaling 
part is negative and at any extremum with negative v the 
scaling part is positive. Since for cr > — 1 the scaling part 
is negative at p = 0 and since the potential (of a stable 
theory) is positive for large values of the field we conclude 
that every globally defined fixed point solution v(p) is 
positive. If v takes negative values for certain fields, then 
it would have a minimum p m \ n with negative v and the 
scaling part would change sign between p = 0 and p m - ln . 
This in turn implies that the scaling part would vanish 
at some field p„ < p m \n or that the potential would be 
singular at p a . 

Critical values of the parameter cr can be detected as 
discontinuities of either p a or its first derivative d a pa • In 
general, these are directly recognizable when plotting the 
function (p a as either jumps or spikes. In the universality 
classes of the scalar model they always appear as jumps 
of p a . Let us illustrate how the critical solutions are de¬ 
tected in practice by solving (15) in the case of d = 3. We 
begin by numerically integrating (15) with initial condi¬ 
tion (19) as a function of the parameter cr, with both 
formulas specialized to the case d = 3 and r) = 0. For 
almost all values of cr, the numerical integration comes 
to an end at a singularity p a defined by ( fl8[ and whose 
value depends on a as shown in Fig. [T| The plot of p a 



FIG. 1. Lo cation of the singularity p a as a function of cr by 
solving (151 for d = 3 and y = 0. 


reveals various salient features. The singularity located 
at the point cr = —1 is due to our choice (19) for the 
parametrization of the initial condition and simply rep¬ 
resents an infinite solution v(p) —> oo. The singularity 
located in the center at cr = 0 is the Gaussian solution 
which trivially extends to all values of the field because 
the potential is identically constant v{p) = const. On 
the right of the plot the function p a slowly increases to 
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infinity where another trivial solution, representing an 
infinite mass limit is located. The interesting feature of 
Fig- [ 1 ] is the singularity located on the left of the Gaus¬ 
sian one at a CI ss —0.1861, which is characterized by a 
marked singular behavior of both (p a (a jump) and its 
first derivative with respect to a (a spike). We can try 
to interpret the meaning of this singularity by observing 
the behavior of the solution in the vicinity of the critical 
value cr cr ~ —0.1861 as shown in Fig. [2j The solutions 



v(<p) 


V 

FIG. 2. Almost critical solutions for the values a = —0.1860 
and a = —0.1862. 


with a slightly on the left of cr cr display a minimum and 
end at the singularity ( [l 8 | with positive values of v((p a ) 
and v'(<Pa). Those with a slightly on the right of a cr end 
with negative values of v{<f a ) and v'{(p a ). The simultane¬ 
ous change of signs in (18) is responsible for the sudden 
increase of (p a while crossing a from left to right as seen 
in Fig. |T| 

The left side of the non-trivial singularity in Fig. 0 d is- 
plays solutions possessing a single non-trivial minimum 
before terminating in the singularity. Therefore these 
are solutions with a non-trivial value of the order pa¬ 
rameter ipo such that v'(ipo) = 0. On the right of the 
singularity the solutions sharply bend to negative values 
and, when following a to zero, they continuously connect 
to the Gaussian solution. It is thus tempting to define 
the critical solution as the one obtained from the limit 
a —> ct cr from the left. The value cr cr has to be determined 
numerically and, within a ten digits precision, takes the 
value 


a cr (LPA) = -0.1860662057. 


A closer look at the numerically determined solution re¬ 
veals the critical properties expected for the potential of 
the Ising-class in three dimension as shown in Fig. [3] The 
critical potential in this figure has been calculated with 
the same precision by integrating the ODE © using 
77 = 0 as an input. We refer to the anomalous dimension 
appearing parametrically in (151 as rj m . The consistency 


of the LPA can be further tested by using the solution 
itself to compute the anomalous dimension through the 
corresponding formula (12). The new result for 77 will be 



.. 

0.0 0.1 0.2 0.3 0.4 


v(<p) 


<P 


FIG. 3. Critical solution obtained for the best determination 
of (rer¬ 


an output of the potential’s determination and in general 
will be different from the input parameter. We thus call 
the new anomalous dimension i ] oui . In the case of the po¬ 
tential of Fig. [3] the minimum is located at ipo = 0.2476 
and the use of |l2[) yields ? 7 ou t = 0.1903, thus outlining 
that, strictly speaking, in the LPA the anomalous dimen¬ 
sion is not determined. 


The consistent determination of 77 required by the 
framework of the LPA’ can be achieved by constructing 
an algorithm that iteratively matches the input param¬ 
eter for the anomalous dimension r]i n with the output 
parameter r] out to a desired precision. The algorithm be¬ 
gins with a (reasonable) guess for 77 ; n and calculates the 
corresponding v aci (<p) which will be different from Fig. 
[3] Then it computes rj ou t associated with the critical po¬ 
tential, and finally iterates by choosing a new value of 
?7i„ closer to 77 out . Further details of the algorithms are 
discussed in Appendix [B] The result of the iteration is 


given with great accuracy in Subsect. IIC 


B. Deformations of the simple scalar model 


Once both a critical value u cr of a and the correspond¬ 
ing solution v*(ip) = v acl (ip) are identified, it is possible 
to compute the critical properties of the solution in the 
form of critical exponents. These are obtained by lin¬ 
earizing in the deviations from the critical solution, which 
means inserting the deformed solution 


v(<p)^v*(<-p) + e6v(<p) ( 20 ) 

into © and linearizing in the infinitesimal expansion 
parameter e. Here ko is a reference scale, and 6 plays 
the role of the critical exponent of the deformation. The 
linear order in e gives an ODE for the fluctuation 5v((p) 
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that in canonical form reads 
i (d + 2) 


6v" —c 


d + 2 — r] 


[1 + v*"(<p)] 2 x 


^(0 — d)8v(ip) + —(d — 2 + rj)p5v' 


( 21 ) 


For later use, it proves convenient to perform a transfor¬ 
mation that redefines the deformation as 


Su(p) = exp ^-a J dyy[l + v*"(y)] 2 ^J 8v(p ), 


d + 2d — 2 + rj 

a = —;-;— - > 0 . 

4cd d + 2- rj 


( 22 ) 


such that the linearized flow equation (211 does not con¬ 


tain the first derivative 5u'. The explicit form of the 
transformed ODE for Su(p) is involved and we found it 
not particularly instructive; thus we do not write it down 
explicitly. 

The second derivative v*"(<p) ' n the linear fluctuation 
equation (211 is an even function and we may assume 


that the deformations 5v{p) are even or odd functions 
which means that they have the same parity as v*(p) 
(parity-even) or the opposite parity (parity-odd). Again, 
due to the linearity of (21) we are free to choose the 


normalization of 8v(ip) at will. This leaves us, modulo 
the global normalization, with two possible choices for 
the boundary conditions of the fluctuations at p = 0 


<5u(0) = 1, <Si/(0) = 0, and (23) 

<Su(0) = 0, &/(0) = l, (24) 


for parity-even and parity-odd deformations, respec¬ 
tively. The same is of course true for the transformed 
deformations 8u(p). The spectrum of the critical expo¬ 
nents 9 of the deformations is quantized when the defor¬ 
mations are required to be polynomially bounded, thus 
adding a further condition at large values of the field. We 
denote the quantized spectrum of exponents as {9+ ,9~}, 
where +/— label the even/odd parity of the deformation 
and the index i £ N orders the exponents starting with 
the fluctuations with the highest positive critical expo¬ 
nent downward. 

For the actual determination of the spectrum it is im¬ 
portant to realize that, in our setup, the solution v*(p) is 
determined numerically and exists within a finite range 
only, \p\ < p act < oo. Hence it is necessary to replace the 
boundary condition at infinity with appropriately chosen 
conditions at <p acr to reproduce the quantized spectrum 
of fluctuations. There are two natural ways to do so. 

On the one hand the requirement that the fluctua¬ 
tions are subleading at infinity can be translated, within 
numerical accuracy, to the requirement the fluctuations 
must tend to zero at the actual numerical boundary 
5v(if aci ) = 0. This can be achieved numerically with the 
help of the shooting method. The errors of introducing 
artificial boundary conditions for 8v(p) can be controlled 
when one relaxes the condition by setting Sv(pa CT ) = c, 


where c is an arbitrary numerical constant. Typically, the 
critical exponents show a very mild dependence (of order 
10 ~ 3 ) on c even on the very large range c £ [—10 3 ,10 3 ]. 
This is due to the fact that the factor [l+v*"(<^)] appear¬ 
ing in (21) mitigates the effects of any boundary condi¬ 
tion by suppressing it in the limit p —> p „ ci . The depen¬ 
dence on the approximate boundary condition can be re¬ 
laxed even further by considering the transformed fluctu¬ 
ation (22), for which it is natural to require Su(p acT ) = 0. 

On the other hand, since the numerical solution exists 
only in a finite range, it seems natural to impose periodic 
boundary conditions for the even fluctuations and an- 
tiperiodic for the odd ones. We studied both cases by im¬ 
plementing the very efficient SLAC derivative method for 
the transformed fluctuations defined in (22) to compute 
the critical exponents. In this case, the potential in the 
range p £ (— p a „, ^a CI ) is discretized over a equidistant 
grid such that the Schrodinger-type differential equation 
for Su(p) turns into a (large) matrix-equation in ip space. 
Most importantly, the derivative with respect to p is dis¬ 
cretized as non-local SLAC derivative. In the past this 
pseudo-spectral method based on trigonometric functions 
proved to be very efficient to accurately calculate the 
low-lying energies of Schrodinger-type operators. One 
may also estimate the characteristic width of an eigen¬ 
mode Su(p), which can be compared with the maximal 
field (p aci to decide whether the solution can be trusted 
and to what extent the spectrum is affected by finite size 
effects. Furthermore, the convergence of the spectrum 
can be checked by increasing the number of grid points. 
Not unexpectedly, the critical exponents obtained with 
the spectral method yield the most accurate results. We 
refer to Appendix [P] for more details on the method. 


C. The Ising class in d = 3 


In this Subsection we provide tabulated results for the 
critical solution of (151 in three dimensions that corre¬ 


sponds to the Ising universality class. In each case we 
studied the critical value cr cr with at least 10 digits pre¬ 
cision, while the anomalous dimensions were determined 
with at least seven digits 

er cr (LPA ) = -0.1860642495; 

o- cr (LPA') =-0.1356610909, 77 = 0.1119795; (25) 

f7 cr (LPA / ) = -0.1657407105, rf = 0.04427234. 

The critical exponents were always determined with at 
least six digits precision. In Table [T] the first three rows 
correspond to the LPA schemes described in the previous 
Sections, including the LPA for which the input anoma¬ 
lous dimension ? 7 ; n = 0, the LPA’ obtained by matching 
?/i n with the anomalous dimension at the minimum using 
(12), and the LPA’ obtained by matching r]i n with the 
anomalous dimension at the minimum using (flS 


The 

fourth row gives the result from a LPA in which the best 
known value of the anomalous dimension is taken from 
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[ 2 j as an input parameter, while the last three rows con¬ 
tain the best known results from various techniques [ 2 ], 
from lattice’s finite-size scaling [M], and from the confor¬ 
mal bootstrap approach [35]. We provide values for the 
critical parameter er cr , for both input and output anoma¬ 
lous dimensions, as well as for the critical exponent v 
governing the scaling of the correlation length which is 
defined as the inverse of the largest even critical expo¬ 
nent v = 1/9+. We also give numerical estimations for 
the dimensionless coupling, defined as 


,( 4 ) 


9cr = 


(¥>o) 




(26) 


Note that the definition (|26j) refers to renormalized quan¬ 
tities and bears an intrinsic local potential scheme de¬ 
pendence. The comparison of this critical coupling with 
analog definitions coming from a lattice’s bare action (see 
for example EH and references therein) might not be 
straightforward. The explicit relation between (26) and 


its lattice’s counterpart could perhaps be computed along 
the lines of a non-perturbative generalization of [35]. All 



O'er 

9 cr 

Vin 

Tjout 

V 

LPA 

-0.1861 

16.1912 

0 

0.1903 

0.6496 

LPA’ ?7 

-0.1357 

9.6098 

0.1120 

0.1120 

0.6453 

LPA’ r/ 

-0.1657 

13.1693 

0.0443 

0.0443 

0.6473 

best r/ 

-0.1693 

13.6776 

0.0364 

0.1626 

0.6477 

p.&v. m 



0.0364 


0.6301 

Lattice [35] 

- 

- 

0.0363 

- 

0.6300 

Bootstrap [36J 



0.0363 


0.6300 


TABLE I. Critical parameters for the Ising class in d = 3. 
The parameters are explained in the main text. The shooting 
and SLAC-method yield the same values for v. 


numbers of Tab. [T] are rounded to 10 4 . 

The non-smooth regulator § might lead to problems 
when it is used beyond the next to leading order of the 
derivative expansion [28) . Thus, one may ascribe part of 
the error in determining the anomalous dimension to this 
particular regulator instead of the LPA’-approximation. 
To investigate this question we repeated the calculation 
with a power law regulator 

K k (q 2 ) = (27) 

Q 

and with the exponential regulators 

W) = ^I. y = fJ ( 28 ) 


which are compatible with the derivative expansion. The 
spike-plots look very similar to those for the non-smooth 
regulator. The LPA’-results with 77 defined in (12) are 
listed in Table El We see that the anomalous dimension 
does not improve considerably compared to the results 
for the non-smooth cutoff pi). It does improve a bit for 
the particular choice (a, b, c) = ( 1 , 1 , 1 ) of the param¬ 
eters in the exponential cutoff. But at the same time 


cutoff 

O'er 

fler 

V 

V 

power law 

-0.3673 

1.5528 

0.1266 

0.6505 

exponentiali 

-0.3068 

9.4133 

0.1197 

0.6462 

exponential 

-0.2808 

12.6742 

0.0690 

0.6530 


TABLE II. The same critical parameters obtained with power 
law and exponential cutoffs. The results in the row marked 
by exponentiali are obtained with parameters (a, b, c) = 
(l,log(2), 1.44) (which are optimal for the LPA) while those 
in the line exponential are obtained with (a, b, c) = (1,1,1). 
The cutoffs are given in the main text. 


the result for the critical exponent v moves away from 
the accurate lattice-result v = 0.6300 quoted in Table |l] 
We conclude that the errors in determining the anoma¬ 
lous dimension must be ascribed mainly to the LPA’- 
approximation. Thus we proceed with the non-smooth 
regulator Q for which the occurring momentum inte¬ 
grals can be calculating explicitly. 

Next we give in Table |Hl| the calculated values for the 
critical exponents corresponding to even deformations of 
the critical Ising solution. For the various local potential 
schemes, we adopted both the shooting method which en¬ 
forces the boundary condition 6v(tp) = 0 and the spectral 
SLAC-method 0 


Shooting: 


0+ 

9 t 


9 t 

LPA 

1.5395 

-0.6557 

-3.1800 

-5.9122 

-8.7961 

LPA’ T) 

1.5496 

-0.5323 

-2.9017 

-5.4604 

-8.1642 

LPA’ r/ 

1.5448 

-0.6076 

-3.0729 

-5.7415 

-8.5789 

best rj 

1.5440 

-0.6162 

-3.0919 

-5.7699 

-8.5993 

SLAC: 

ot 

6+ 

9 t 


9 i 

LPA 

1.5395 

-0.6557 

-3.1800 

-5.9122 

-8.7961 

LPA’ rj 

1.5496 

-0.5323 

-2.9017 

-5.4598 

-8.1585 

LPA’ r/ 

1.5448 

-0.6076 

-3.0727 

-5.7384 

-8.5513 

best rj 

1.5440 

-0.6163 

-3.0920 

-5.7698 

-8.5956 


TABLE III. Critical exponents of the even deformations for 
the Ising class in d = 3. 


Finally, we give in Table |IV| the numerical values for 
the critical exponents corresponding to odd deformations 
of the Ising scaling solution, again displaying the results 
of both shooting and spectral SLAC-method. 


2 Here and for the rest of the paper we do not display the critical 
exponent ()/ corresponding to a constant deformation of the po¬ 
tential that determines the scaling of the volume operator, since 
for all critical models it is equal to the dimensionality of the 

S stem (which in this example is three). We refer to Appendix 
for more details on 0/ and its peculiar behavior away from 
criticality. 
















































D. The Ising class in d = 2 


Shooting: 

01 

O2 

O3 

01 

O5 

LPA 

2.5000 

0.5000 

-1.8867 

-4.5244 

-7.3376 

LPA’ tj 

2.4440 

0.5560 

-1.6892 

-4.1635 

-6.8188 

LPA’ 7 / 

2.4779 

0.5221 

-1.8103 

-4.3854 

-7.1387 

best 77 

2.4818 

0.5182 

-1.8240 

-4.4097 

-7.1665 

SLAC: 

0i 

O2 

63 

01 

O5 

LPA 

2.5000 

0.5000 

-1.8867 

-4.5244 

-7.3377 

LPA’ 77 

2.4440 

0.5560 

-1.6891 

-4.1608 

-6.7933 

LPA’ 77 ' 

2.4779 

0.5221 

-1.8102 

-4.3845 

-7.1288 

best 77 

2.4818 

0.5182 

-1.8240 

-4.4097 

-7.1665 


TABLE IV. Critical exponents of the odd deformations for 
the Ising class in d = 3. 


As expected, in all cases we detect one relevant defor¬ 
mation having the same even parity as the critical so¬ 
lution. Two more relevant deformations show up in the 
odd sector and can be related to the anomalous dimen¬ 
sion ? 7 i n which governs the scaling part as explained in 
Appendix [A] In particular, the first two exponents of the 
odd sector have the numerical values 9^ ~ (5 — 77 i n )/2 
and O 2 — (1 + 77 i n )/2 consistently with the results of Ap¬ 
pendix [A] 


Our estimates for the critical exponent v deviate at 
most by 3 percent from the most accurate values com¬ 
piled in [2]. The agreement is, however, not as excel¬ 
lent for the anomalous dimension and, in particular, for 
the LPA’ scheme based on the definition (12 1 , which dif¬ 
fers from the best known value by a factor three. This 
discrepancy is a known limitation of the local potential 
schemes and could be eased by improving the LPA to 
a higher order derivative expansion [22H24] . We discuss 
this point further in Sect. [Vj 


For completeness we also list the critical exponents 
for the even and odd deformations as obtained with the 
power law regulator (271 and the exponential regulator 
(28) in Table [V] The parameters of the exponential reg¬ 
ulator are the same as in Table |TT| and again the anoma¬ 
lous dimension is determined with the formula (12). The 
exponents are consistent with those calculated with the 


non-smooth regulators and listed in Tables III and IV 


cutoff 

oi 

O2 

of 

oi 

oi 

power law 

1.5373 

-0.4954 

-2.7695 

-5.1950 

-7.7301 

exponentiali 

1.5475 

-0.5207 

-2.8689 

-5.4011 

-8.0726 

exponential 

1.5313 

-0.5671 

-2.9482 

-5.5109 

-8.2175 

cutoff 


^2 

O3 

O3 

O5 

power law 

2.4367 

0.5633 

-1.6095 

-3.9665 

-6.4505 

exponentiali 

2.4402 

0.5598 

-1.6676 

-4.1155 

-6.7210 

exponential 

2.4655 

0.5345 

-1.7301 

-4.2101 

-6.8468 


TABLE V. Critical exponents of the even and odd deforma¬ 
tions for the Ising class in d = 3 and calculated with a power- 
law and two exponential regulators. The parameters are the 
same as in Table El 


In two dimensions and for zero anomalous dimension 
the fixed point equation © can be integrated without 
ever encountering a singularity for all initial conditions. 
The resulting solutions have been argued to belong to 
the Sine-Gordon universality class [16], meaning that the 
LPA is not capable of detecting the two-dimensional crit¬ 
ical solutions known as minimal models. To circumvent 
this deficiency one may continue (15) to a dimensional¬ 
ity slightly above two d > 2 , or allow for a non-trivial 
anomalous dimension rji n . In both cases, the effective di¬ 
mensionality of the held (d — 2 + rj ln )/2 : which vanished 
in LPA, becomes positive. Fixing d = 2 and using a non¬ 
zero 77 i n , it is possible to compute the function (p a as in 
the previous section. With decreasing r\- m the correspond¬ 
ing plot of <t —^ (p a shows an increasing number of singu¬ 
larities, which we identify with the well-known minimal 
conformal models in two dimensions. Fig. [4] shows the 
spike-plot of Cp a for rj[ n = 1/4. It reveals that four critical 
solutions are already present for this choice of r). In the 
limit of effective field’s dimension (d— 2 + r/i n )/2 —► 0 , so¬ 
lutions appear spawning from the Gaussian one at a = 0 
and alternating from left to right. A closer inspection re- 



FIG. 4. Location of the singularities (p a as a function of a 
by solving ( |15| ) for d = 2 and r/ = 1/4. 

veals that the leftmost non-trivial singularity corresponds 
to an Ising class, the rightmost to a tri-Ising class, the sec¬ 
ond from the left to a quadri-Ising and the second from 
the right to a penta-Ising. It appears that the heights 
of the singularities, that is their approximated domain 
of existence, are ordered like the inverse of the central 
charge of the corresponding minimal model. 

We now concentrate our attention on the LPA’ of the 
two-dimensional Ising class; that is, we isolate the left¬ 
most singularity of Fig.[4]and match the input and output 
anomalous dimensions. Like in the previous Subsections, 
this can be done by either using 77 or the Goldstone’s 77 ', 
resulting in two different sets of critical parameters 

cr cr (LPA') = -0.3583283424, 77 = 0.436348; 

cr cr (LPA') = -0.5291046049, ?/ = 0.233594. 


( 29 ) 
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Analogously to (26) we define the two-dimensional criti¬ 
cal coupling as the ratio 


~~ 6v"((p 0 ) ’ 


(30) 


which corresponds to the canonically dimensionless ratio 
g cr = A/to 2 of the quartic coupling and squared mass in 
two dimensions. 

In Tab. |VI| we listed the critical parameters for the 
Ising class in two dimensions. The first two rows contain 
the LPA’-values obtained with the anomalous dimensions 
?7 and rf. In the third row g- ln is set equal to the exact 
value taken from CFT, and finally, the last row contains 
the exact results coming from CFT. 



O'er 

9cr 

7?in 

9out 

^shoot 

^slac 

LPA’ g 

-0.3583 

13.0397 

0.4363 

0.4363 

1.0550 

1.0551 

LPA’ 77 ' 

-0.5291 

55.5004 

0.2336 

0.2336 

1.3265 

1.3264 

Exact 77 

-0.5143 

47.9118 

0.2500 

0.7477 

1.2871 

1.2870 

CFT 



0.2500 


1.0000 

1.0000 


provides the least accurate value for the anomalous di¬ 
mension. We refer to the discussion in Sect. 0 for more 
details concerning this point. 


E. The TVi-Ising class in d = 2 


In this Subsection we briefly summarize our findings 
about the rightmost singularity appearing in the spike- 
plot associated to the tricritical Ising class in Fig. [3] The 
properties of this class can be investigated by adopting 
the same techniques that were used for the Ising class. 
Again, the critical parameters can be determined by ei¬ 
ther using g or the Goldstone’s rf for the matching of g- ln 
and r] out , thus resulting in two different sets of critical 
parameters 


cr cr (LPA') = 0.2597409360, 77 = 0.311938 ; 
cr cr (LPA') = 2.9267843351, 77 ' = 0.0686649. 


(31) 


We also supplement the critical coupling g CI in (301 by 


TABLE VI. Critical parameters for the Ising class in d = 2. 


In Tab.s [Vn| and |VTTT| we listed the critical exponents 
for parity even and odd deformations of the critical so¬ 
lutions. As before, we provide both the results obtained 
with the shooting and the spectral SLAC-methods. 


Shooting: 

si 

St 

St 

si 

si 

LPA’ 77 

0.9478 

-0.6761 

-2.6194 

- 4.7373 

-6.9765 

LPA’ g' 

0.7539 

-0.8935 

-2.9866 

-5.3007 

-7.7637 

Exact g 

0.7770 

-0.8741 

-2.9597 

-5.2613 

-7.7111 

SLAC: 

si 

st 

St 

si 

si 

LPA’ g 

0.9478 

-0.6761 

-2.6193 

-4.7347 

-6.9938 

LPA’ g' 

0.7539 

-0.8935 

-2.9863 

-5.2960 

-7.7125 

Exact g 

0.7770 

-0.8741 

-2.9593 

-5.2547 

-7.6424 


TABLE VII. Critical exponents of the even deformations for 
the Ising class in d = 2. 


Shooting: 

Si 

Si 

S3 

Si 

S5 

LPA’ g 

1.7818 

0.2182 

-1.6189 

-3.6605 

-5.8438 

LPA’ g' 

1.8832 

0.1168 

-1.8903 

-4.1204 

-6.5150 

Exact g 

1.8750 

0.1250 

-1.8703 

-4.0880 

-6.4689 

SLAC: 

Si 

Si 

S3 

Si 

S5 

LPA’ g 

1.7818 

0.2182 

-1.6189 

-3.6606 

-5.8452 

LPA’ g' 

1.8832 

0.1168 

-1.8903 

-4.1206 

-6.5168 

Exact g 

1.8750 

0.1250 

-1.8703 

-4.0880 

-6.4696 


TABLE VIII. Critical exponents of the odd deformations for 
the Ising class in d = 2. 


The estimates for the critical exponent of the corre¬ 
lation length in Tab. |VI| deviate from the exact values 
determined within CFT by 6 % — 30 % with the standard 
LPA’ scheme giving the best agreement. On the other 
hand and similarly to Subsetc. IIC the LPA’ scheme 


= a (6) (y>o) 

56,cr “ 120 v"(tpo) 


(32) 


which corresponds to the dimensionless ratio <? 6 ,cr = 
Ae/ to 2 for the expansion 


1 \ mr 2 , A 4 Ag 6 

+ ~^P + +■■■ 


(33) 


accommodating a higher criticality for the tricritical Ising 
class. 

The critical parameters are determined within the var¬ 
ious local potential schemes in Tab. m In Tab.s [X] and 
XI we give the critical exponents for parity even and odd 
deformations. We provide only the results obtained with 
the spectral SLAC-method, but we also compare with 
the exact results coming from CFT. 



O cr 

9cr 

P 6 ,cr 

gin 

9out 

V 

LPA’ g 

0.2597 

8.4410 

23.1728 

0.3119 

0.3119 

0.6001 

LPA’ g' 

2.9268 

182.135 

19210 

0.0687 

0.0687 

0.5239 

Exact g 

1.0639 

38.8756 

747.039 

0.1500 

0.7256 

0.5498 

CFT 

- 

- 

- 

0.1500 

- 

0.5560 


TABLE IX. Critical parameters for the tricritical Ising class 
in d = 2. 


SLAC: 

si 

St 

si 

St 

St 

LPA’ g 

1.6665 

0.6836 

-0.4510 

-1.7920 

-3.2840 

LPA’ g' 

1.9087 

0.3596 

-0.7899 

-2.1364 

-3.7899 

Exact g 

1.8187 

0.5501 

-0.6248 

-2.0554 

-3.6878 


TABLE X. Critical exponents of the even deformations for 
the tricritical Ising class in d = 2. 


It is worth noting that our definition of the critical cou¬ 
plings (321 displays a wide range of values. This is espe¬ 
cially true for the quantity ge >CI which varies over several 
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SLAG: 

*r 

^2 

^3 

81 

8s 

LPA’ g 

1.8440 

1.1560 

0.1560 

-1.1008 

-2.5218 

LPA’ g' 

1.9657 

0.6890 

0.03434 

-0.7899 

-2.1364 

Exact g 

1.9250 

0.9715 

0.07500 

-1.3157 

-2.8562 


TABLE XI. Critical exponents of the odd deformations for 
the tricritical Ising class in d = 2. 


orders of magnitude between the various LPA schemes. 
The numerical values of the parameters g cr and g§ t cr can 
be used in (331 to obtain an approximate form of the 
local potential solutions. Then a decreasing anomalous 
dimension in Tab. m leads to a solution with a more 
pronounced minimum, which means that the coefficients 
in the Taylor-expansion have large numerical values. 


III. THE JV = 1 WESS-ZUMINO MODEL 

The most straightforward supersymmetric extension 
of the scalar field theory is the supersymmetric JV = 1 
Wess-Zumino model (WZM). Since in 3-dimensional Eu¬ 
clidean space there are no Majorana-spinors we formulate 
the flow equations in Minkowski space. Then Majorana- 
spinors exist in three and two dimensions and have two 
components in both Minkowski spaces. Hence an interpo¬ 
lation of the flow equations from three to two dimensions 
should be possible. Both in three and two dimensions the 
field content of the model can be characterized through 
the superfield 


$ = <p+(tol>) + ±(00)F, (34) 

in which p is a real scalar field, V a Majorana spinor 
and F a real auxiliary field. We defined the Grassmann 
variable 9 and its conjugate 9 which will play the role of 
coordinates in superspace. 

Following previous works [50H52] we formulate the LPA 
by using a manifestly supersymmetric truncation for the 
effective average action in superspace, 


on-shell effective average action 

ponsM! = f d d X (^ d ^p + W FM 

J v ^ ( 37 ) 

- - iw'( , (‘p)('0V’)) • 

A direct comparison of the first line with ([5]) shows that, 
on-shell, the scalar sector of the WZM is equivalent to 
a scalar field theory controlled by a potential V k (p) = 
W' 2 {p)/2Z k and coupled to the Majorana field with a 
particular Yukawa interaction whose form is constrained 
by supersymmetry. 

It is possible to construct an infrared cutoff like © 
which is quadratic in the field multiplet {p^ip^F} and 
that preserves supersymmetry, thus giving rise to a renor¬ 
malization group flow that is fully covariant with respect 
to supersymmetry. We refer to for more details 

on the construction of a manifest supersymmetric flow, 
while we give the simplest implementation that (after 
Wick-rotation to Euclidean space) in momentum space 
is 


Zu 

AS k = yr 
2 


J d d q(p q TZ k (q 2 )p- q + F c 


K k (q 2 ) 


F_ n 


1 7 ^ K k (q 2 ) 


'I’-q) > 


(38) 


which depends on a single cutoff kernel lZ k (q 2 ) carry¬ 
ing the scheme dependence of the flow. Notice that 
a manifestly supersymmetric RG-flow with a regulator 
quadratic in the fields can only be formulated in the off- 
shell formulation. The auxiliary field can only be elimi¬ 
nated after the complete FRG-flow has been determined. 
Nonetheless, we still refer to V k (p) = W k (p)/2Z k as the 
scale-dependent effective potential of the WZ-model, re¬ 
membering that in the limit k —> 0 it does coincide with 
the effective potential for the scalar mode p. 

As for the simple scalar model it is possible to explic¬ 
itly evaluate the RG flow of the model within the given 
truncation once the cutoff kernel 7 Z k (q 2 ) is explicitly cho¬ 
sen. A simple family of kernels that mimics © is 


T k = j d d x J d9dd(~Z k <f>K$ + 2 Wfc($)) , (35) 


in which we introduced the second order superderivative 
operator I\. This truncation generalizes ([5]) by introduc¬ 
ing the scale-dependent effective superpotential W k {<&). 
Integrating explicitly over the Grassmann variables ( |35[ ) 
can be expressed in terms of the component fields 


r k 


d d x(^^d tl pd^p - l -Z k V$V + -^ p2 
- l -WM{^) + W' k {p)F). 


(36) 


It is instructive to eliminate the auxiliary field F through 
its algebraic equation of motion to obtain the so-called 


F k {q 2 ) = q 2 ( ( k / \q\) n - l) 9(k 2 - q 2 ) , (39) 

which is parametrized by the exponent n. The choice 
n = 2 leads to a finite flow both in two and three di¬ 
mensions, while the choice n = 1, which coincides with 
©, yields a divergent flow in two dimensions. In com¬ 
plete analogy to the scalar field theory, we introduce the 
dimensionless renormalized field = Z^^k^-^^p 
that is the natural argument of the dimensionless renor¬ 
malized superpotential 

w k (p R ) = k~ d / 2 W k (p)/Z k , (40) 

with anomalous dimension g = —kd k Z k /Z k . Once again, 
for notational simplicity we omit the labels in the fol¬ 
lowing. In the evaluation of the Wick-rotated flow Q 
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within the truncation (35) the dependence on the space- 
dimension is easily seen in momentum space. This way 
one arrives at a d-dependent flow for the superpotential 
which serves as a natural analytic continuation of the RG 
flow of the model to dimensions 2 < d < 3. 

The flow is computed in practice by projecting (W) to 
the simple field configuration <f> = tp = const. Using ( 39] ) 
for n = 2 we obtain the system 


kd k w(ip) = S w>r ,[w,w';(p] + T w ^[w "), 
S Wt r,[w,w'-,ip] = (1 - d)w(ip) + -—, 


(41) 


di] 

(d + 2) 

where we maintained a compact notation by defining 


= -y ((2 - ^i,oM + 


w 


Ai tm [w"] = 2 Fi (l, m + d/2, m + 1 + d/2, —w" 2 ) 


= \m + 


r» 1 j.m+d/2 —1 

d t 


(42) 


(1 + tw" 2 ) 1 


which depends on the Gauss hypergeometric function 
2 T|. The integral representation shows explicitly that 
for non-negative l and m the function Ai >m ]w"\ is a pos¬ 
itive function on the real axis which satisfies 


Al, m [-w"] = Al, m [w"], Al,m[ 0] = 1 . (43) 

It follows that J~ w . r ) is an odd function of w" which is 
negative for positive arguments. Its negative slope at the 
origin is 


•T, [ o] = -<* (i - jh) ■ (44) 

In an off-shell formulation of supersymmetric field the¬ 
ories there is an efficient way to extract the anomalous 
dimension, see [M1I3T1 39] for details. For consistency the 
anomalous dimension must be evaluated at a particular 
field configuration, which we choose to be the order pa¬ 
rameter <po. Introducing the dimensionless renormalized 
on-shell potential 


__ w’(ip) 2 

v{<p) = -O- 


(45) 


the order parameter is defined as its (nontrivial) mini¬ 
mum through the relation v'(<pa) = 0. Using the com¬ 
pact notations Wq = w"(<po) and w= w"'(<po), the 
anomalous dimension can be expressed implicitly as 


V = Cd 


(2 - 77)^43,0^0] + ( f+2^ 3 ’ 1 [ w o] 

^(2 77 ) . r //2 /// 

--43,lKJ w 0 w 0 


d + 2 
dr] 


(46) 


d + 4 


a r //1 //^ 

A 3 , 2 [W 0 \W 0 


In the limit d j. 2 the equations (|41|) and (46) can be seen 
to reproduce the results in [301 131|. Analogous results 


are obtained for the choice n = 1 in (39) that are finite 


above 2 dimensions, and in the limit d f 3 reproduce the 
results of [32] - Our computation with the cutoff n = 2 
offers a consistent analytic continuation of the flow of 
the superpotential which interpolates between the known 
results for d = 2,3. The results obtained with the n = 1 
cutoff will only be used in 3 dimensions to determine, to 
some extent, the cutoff independence of the methods. 

To conclude this Section, let us note that, in contrast 
to the non-supersymmetric case, the anomalous dimen¬ 
sion (46) has, for some models, a nonzero value when 


evaluated at the configuration <po = 0. Let us thus define 
a further anomalous dimension 


Vo = Vl^o 


(47) 


For all the critical solutions that we are interested in, the 
location ipo = 0 is either a minimum or a local maximum, 
but never the global minimum which is needed in the 
LPA’. This anomalous dimension is thus perturbative, in 
the sense that it is not obtained by an expansion about 
the true vacuum state. Nevertheless, for the WZ-model 
it is possible to investigate a new local potential scheme, 
which we call LPA°, that uses rjo as an input for the 
integration of (48). This partly perturbative scheme lies 
in between the LPA and the LPA’. 0 


A. Scaling solutions for the WZM 


It seems difficult to explicitly rearrange the ODE cor¬ 
responding to the stationary solutions of the highly non¬ 
linear flow equation (41) in a standard form w" = ... as 


we did for the scalar field in (15). Rather than finding 


stationary solutions for w(tp) via ( |4l| ), we decided to in¬ 
vestigate the stationary solutions of the flow equation for 
w'(ip), which can be cast in standard formj^] To obtain 
this equation we differentiate (41) with respect to tp and 


solve the fixed-point equation kd k w'(p) = 0 for w"'(ip) 
with the result 


'(¥>) = 


1 


2 


d - 77 , d- 2 
w -— 


V 


-tpw 


')■ (48) 


However, the asymptotic structure of stationary solutions 
is more directly (and more comfortably) accessible from 
(41) itself. For large values of w" the function T w ^ ap¬ 
pearing on the right hand side of ([47]) has the asymptotic 


3 In passing we note that this scheme was originally dubbed as 
LPA’ in 13011311 , where a different definition for the local potential 
approximation for the WZ-model has been used. 

4 It is actually possible to numerically solve © by applying an 
implicit solver combined with a residual method such as New¬ 
ton’s. This however comes at the expense of the solution’s accu¬ 
racy. We explicitly checked that the results reported below can 
be reproduced by the sole use of GD- 
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expansion 


J w,rj — 


^+( 2 - 7 ?)—^- 
2w" V ''d-2 

C2 


2 < d < 4; 


(77 + (2-7;)log(u/ /2 )), d = 2, 


(49) 

which implies, like in the scalar field theory, a divergent 
second derivative of the superpotential when the scaling 
part of the equation is zero. Neglecting the leading logs 
and overall factors 




1 


S„ 


, when S wr >~0. 


(50) 


A closer inspection of (481 reveals that the ODE might 


possess a singularity that is met whenever T' w rj = 0. This 
is indeed true because the function T w ^ r] intersects the 
u>"-axis at the origin with negative slope, see (44). Since 
in addition J r WJI ['w"] tends to zero for large arguments, 


see (49), we conclude that there must exist a minimum 


w 


of T w , v at which T' wrj vanishes. When integrating 
(48) we expect to find a singular w'" when w" approaches 
* In other words, for generic initial conditions we find 


w 


w’"(ip)~ oo, when IF' wr] {w") —> 0 , (51) 

therefore there exists a singular value of the field (p for 
which 


w"(<p) 


when tp 


tp. 


(52) 


A comparison of the definition (521 with the analog def¬ 


inition for the simple scalar model (18) reveals that the 


integrations of the stationary ODEs meet the singulari¬ 
ties in two rather different ways. It turns out that the 
singularity (52) can in fact be crossed if 


d _ d-rj . d-2 
dv Sw ’ v ~ 2 W + 2 




ipw" = 0 


(53) 


at <p = ip, which will turn out to be true only for all 
the critical solutions that we will investigate. When this 
condition is met at a critical point, the third derivative 
of the superpotential does not diverge according to ( [FT] ) 
and the solution acquires the desired global nature. For 
more details on this extension we refer to Appendix [EJ 
The behavior of the singularities must be studied as a 
function of the initial conditions of the fixed point equa¬ 


tion (481. Our only constraint is that the on-shell poten¬ 
tial must be Z 2 -invariant. Recalling that T' w v in (48) is 


an even function of w" it is easy to see that the function 
w'(tp) must either be even or odd in tp. Ifw'(ip ) is even 
then w"(0) = 0 and we parametrize the initial condition 
with w'"(0) = (- so that 


«'( 0 ) = - 

V'(o) =o, 


2 c d (d - rj + 2) 
(d + 2 )(d — rj) 


C-, 


(54) 


which thus serve as boundary conditions for parity odd 
superpotentials. If instead w r (ip) is odd then x</(0) = 0 
and we choose w" (0) = £ + so that 


i/(0) =0, 

«"(0) =C+, 


(55) 


which are the appropriate boundary conditions for parity 
even superpotentials. 

The situation is different from the previously consid¬ 
ered scalar held theory in which one parameter a in 
(19) was sufficient to parametrize all possible critical 
models. In the supersymmetric case the solutions fall 
into two distinct classes characterized by different initial 
conditions (54) and (55) with corresponding parameters 
> 0. For the ease of numerical comparison between 
the scalar-model LPA and the WZ-model LPA, we can 
however provide, using the on-shell dimensionless poten¬ 
tial v(ip) = w' 2 (ip)/2 , the value of cr = n"(0) which is 
defined in analogy to the previous Sections. We obtain 



2c d (d-r] + 2) 
(•d + 2 )(d - r)) 


C_ <0, 


or 


(56) 


= C+ > 0. 


for even or odd superpotentials respectively. This shows 
that the first class will contain solutions with an even 
number of minima as it happens for cr < 0, while the 
second class contains solutions with an odd number of 
minima. We dedicate the next Subsections to the critical 
models that can be detected by plotting the maximum 
value of the field (52) as a function of £±. 


B. The SUSY-Ising class in d = 3 

The supersymmetric extension of the Ising class in 3 di¬ 
mensions is obtained by integrating the fixed point equa¬ 


tion (48) with the boundary conditions (54) for d = 3. 


It is sufficient to study the range of initial conditions for 
which £_ > 0, since the negative can be reached by 
mirror symmetry. Within this range, all numerical solu¬ 
tions of the fixed point equation end at a singularity ipg_ 
obeying ( [52]) such that w" diverges. The plot of (p^_ is 
given in Fig. [5] and should remind the reader of the nega¬ 
tive a range of Fig. [I] as in fact the parameters a and 
can be explicitly related through (56). Even though the 
jump at the critical value C-,cr is not as pronounced as 
in the non-supersymmetric case, the plot of Fig. [5] clearly 
displays a singularity at which the maximal field ip^_ is 
discontinuous and its first derivative diverges. The sin¬ 
gularity is not always as pronounced, and in particular 
it is absent in the LPA scheme with the n = 2 cutoff. 
We therefore decided to display the plot using the LPA 0 
scheme defined by (47), being the closest analogue to the 
LPA. This, in particular, implies that, by varying in 
Fig. [5] both the value of w'( 0) and of ??o vary. 

A closer inspection of the solutions in the vicinity of 
the singularity reveals that the solutions with larger 
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FIG. 5. Location of the singularity of as a function of 
£_ by solving ( |48| ) (LPAo) for d = 3 and g = go using the 
n — 1 cutoff. 



C— jCr 

O cr 

9cr 

Pin 

'Tjout 

V 

LPA u n =i 

2.0239 

-0.0692 

13.5783 

0.1880 

0.1599 

0.7112 

LPA’„=i 

2.1333 

-0.0384 

10.4531 

0.1736 

0.1736 

0.7076 

LPA u n =2 

2.3558 

-0.0641 

10.1464 

0.1807 

0.1531 

0.7094 

LPA’„=2 

2.4794 

-0.0708 

11.0943 

0.1670 

0.1670 

0.7060 


TABLE XII. Critical parameters for the supersymmetric Ising 
class in d = 3. 


than the critical value have an Ising-like potential, while 
those with smaller £_ correspond to solutions which bend 
down and, with decreasing value of £-,cr, an additional 
minimum develops. This is in complete analogy with the 
situation described in Subsect. m for the Ising-class, 
especially if we recall that a oc — Cff. We thus investigate 
the supersymmetric Ising solution as the limit C- I C—,cr 
from above. 

Similarly as for the scalar field theory, we now apply 
the numerical algorithms discussed in the previous sec¬ 
tion to investigate the critical properties of the supersym¬ 
metric fixed point solutions. In Table |Xll| we provide the 
critical value for £_, the corresponding critical value a CI 
(using (56)), the value of the critical coupling (which is 
defined analogously to (261 by using the on-shell poten¬ 


tial), the anomalous dimension and the critical exponent 
of the correlation length. These numerical estimates are 
given for the choices n = 1, 2 of the cutoff (391 and for the 
available schemes of LPA° and LPA’. The critical values 
of the parameters are: 


C-(LPA 0 „ =1 ) = 2.023871922, g = 0.188003; 
C-(LPA'„ = i) = 2.133306022 , 77 = 0.173624 ; 
£_(LPA°„ =2 ) = 2.355854024, 77 = 0.180673; 
C_(LPA'„ =2 ) = 2.479354615, 77 = 0.167018. 


(57) 


In Tables XIII and XIV we provide our best estimates 
for the critical exponents of the parity-odd and parity- 
even deformations of the superpotential, which corre¬ 
spond to deformations of the on-shell potential of op¬ 


posite parity. 


Shooting: 

01 

Of 

of 

Of 

of 

LPA u „=i 

1.4060 

-0.3510 

-2.5715 

-5.1730 

-8.1200 

LPA’ n =i 

1.4132 

-0.3824 

-2.6813 

-5.4004 

-8.5042 

lpa u „ =2 

1.4097 

-0.3500 

-2.5281 

-5.0357 

-7.8328 

LPA ’„=2 

1.4165 

-0.3773 

-2.6200 

-5.2216 

-8.1418 

SLAC: 

of 

Of 

of 

Of 

of 

LPA U „ =1 

1.4060 

-0.3510 

-2.5715 

-5.1730 

-8.1200 

LPA’ n =x 

1.4130 

-0.3824 

-2.6813 

-5.4004 

-8.5043 

lpa u „ =2 

1.4097 

-0.3500 

-2.5281 

-5.0357 

-7.8328 

LPA ’„=2 

1.4165 

-0.3773 

-2.6200 

-5.2216 

-8.1418 


TABLE XIII. Critical exponents of the even deformations for 
the on-shell potential of the supersymmetric Ising class in 


d = 3. 


Shooting: 

Of 

Of 

Of 

Of 

Of 

LPA u n =i 

0.5940 

-1.4102 

-3.8274 

-6.6048 

-9.7168 

LPA’„=i 

0.5868 

-1.4759 

-3.9911 

-6.9056 

-10.1945 

LPA u n =2 

0.5903 

-1.3941 

-3.7438 

-6.3997 

-9.3328 

LPA’ n =2 

0.5835 

-1.4500 

-3.8781 

-6.6435 

-9.7144 

SLAC: 

Of 

Of 

Of 

Of 

Of 

LPA u n =i 

0.5940 

-1.4102 

-3.8274 

-6.6048 

-9.7168 

LPA’„=i 

0.5868 

-1.4760 

-3.9916 

-6.9056 

-10.1945 

LPA u n=2 

0.5903 

-1.3941 

-3.7438 

-6.3997 

-9.3328 

LPA’„=2 

0.5835 

-1.4500 

-3.8791 

-6.6435 

-9.7144 


TABLE XIV. Critical exponents of the odd deformations for 
the on-shell potential of the supersymmetric Ising class in 


d = 3. 


The first exponent of Tab. |XIlT| corresponds to the de¬ 
formation Sw(p) ex (p, which, according to the discussion 
given in Appendix [A[ has critical exponent (d — g)/2 , and 
is also related to the critical exponent of the correlation 
length by the relation v = 1/9+. We can thus test that 
the relation 


9 


+ 

1 


1 

v 


d — 77 
2 


(58) 


known as superscaling relation [39] . is satisfied within 
our accuracy. The first critical exponent of Tab. |XIV| 
corresponds instead to the deformation 8wfp) ~ 
which is also discussed in Appendix [A] and whose critical 
exponent is related to the anomalous dimension by the 
formula 


__ d-2 + 77 
1 “ 2 


(59) 


This relation is perfectly satisfied within our accuracy. 


C. The SUSY-Ising class in d = 2 


A supersymmetric extension of the Ising class in 2 di¬ 
mensions is obtained by integrating the fixed point equa¬ 
tion (48) with the boundary conditions (54) for d = 2. 
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The plot of the endpoints of integrations is given within 
the LPA° scheme in Fig. [h] and shows that, as expected, 
many critical points emerge in the two dimensional limit. 
A direct inspection of the solutions and the relation (56) 
reveal that these singularities correspond to even criti¬ 
cal solutions of higher criticality. Apparently, the plot 


Shooting: 

8 J 

8 t 

83 

n 

ef 

LPA u n =a 

0.7807 

-0.4385 

-2.2244 

-4.5811 

-7.5113 

LPA’ n =2 

0.8014 

-0.5132 

-2.5460 

-5.3257 

-8.8585 

SLAC: 

8 ]t 

e+ 


n 

ot 

LPA u b=2 

0.7807 

-0.4383 

-2.2224 

-4.5711 

-7.4777 

LPA’ n =2 

0.8014 

-0.5125 

-2.5416 

-5.3095 

-8.8144 



i- 


TABLE XVI. Critical exponents of the even deformations for 
the on-shell potential of the supersymmetric Ising class in 
d = 2 . 


Shooting: 


82 

83 

84 

81 

LPA U „=2 

0.2193 

-1.2610 

-3.3316 

-5.9739 

-9.1951 

LPA’„=2 

0.1984 

-1.4384 

-3.8423 

-6.9975 

-10.9101 

SLAC: 


82 

83 


85 

lpa u „ =2 

0.2193 

-1.2602 

-3.3268 

- 5.9549 

-9.1393 

LPAA=2 

0.1986 

-1.4365 

-3.8335 

-6.9700 

-10.8432 


TABLE XVII. Critical exponents of the odd deformations for 
the on-shell potential of the supersymmetric Ising class in 
d = 2 . 


FIG. 6. Location of the singularities of (p^_ as a function 
of by solving (481 for d = 2 and r/ = r/o (LPAo) using the 
n = 1 cutoff. The rightmost singularity corresponds to the 
Ising model, and then towards the origin all even solutions 
with higher criticalities appear. 


in Fig. [ 6 ] looks qualitatively different from our previous 
plots such as in Fig. [5j It does not show spike-like sin¬ 
gularities, but rather what appear to be bends or flexes. 
It can, however, be checked that at each singular point 
the function is indeed discontinues with a diverging 
first derivative. The only difference from all the previous 
plots is that the first derivative of (p^_ has the same sign 
on both sides of the singularity. 

Approaching the rightmost singularity from the right 
we obtain a supersymmetric extension of the Ising class. 
The critical values for the two available schemes are 

C-(LPA° n=2 ) = 1.7593594599, 77 = 0.438619; , N 
, (60) 

C- (LPA n=2 ) = 1.9477913122 , 77 = 0.397108. 


The numbers for all critical parameters are listed in Tab. 
|XV| As before we also provide the numbers for the criti¬ 
cal exponents of even and odd deformations of the critical 
solution in Tab.s XVI and XVII| respectively. The first 



C — ,cr 

O' cr 

9 cr 

7?in 

T]out 

V 

lpa u „ =2 

1.7594 

-0.2809 

9.0576 

0.4386 

0.3386 

1.2809 

LPA’ n =2 

1.9478 

-0.3393 

11.6385 

0.3971 

0.3971 

1.2478 


D. The SUSY-TVi-Ising class in d = 2 


We now turn our attention to the boundary condition 


(55) corresponding to even solutions w(ip) of the fixed 
point equation. For this condition it is easy to see that 
770 in (47) vanishes such that LPA 0 becomes LPA. Fur¬ 


thermore, the integration of (48 1 for d = 2 and 77 = 0 


does not meet any singularities in analogy to what hap¬ 
pened in both Sect .s [ITD] and |TirCl because the canonical 
dimension of the field is zero. In Fig. [7] we plot the end¬ 
point of the integration of (48) for d = 2, boundary con¬ 


ditions (55) and trial anomalous dimension 77 = 1/10 as 
a function of the parameter £ + . The plot is qualitatively 



<r + 


TABLE XV. Critical parameters for the supersymmetric Ising 
class in d = 2. 


critical exponents of both Tab.s |XVI| and |XVII| are re¬ 
lated to the anomalous dimension through the relations 
(58) and (59) specialized to d = 2. Both are verified 
well within our accuracy, giving yet another check of the 
superscaling relation. 


FIG. 7. Location of the singularities of (p^ + as a function of 
£+ by solving ( |48| ) for d = 2 and r/ = 1/10 using the n = 1 
cutoff. The rightmost singularity corresponds to a Gaussian 
solution. The next singularity on the left is the tricritical Ising 
model. All other solutions with increasing odd criticality are 
arranged from there towards the origin. 

similar to Fig. |6j displaying the same type of singular- 
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ities in the form of discontinuities of tpc+ across which 
the first derivative does not change sign. A closer inspec¬ 
tion of the solutions shows that the rightmost singularity 
seems to represent a transition between a free theory with 
quadratic potential for large values of £+ and a solution 
with two maxima away from tp = 0 . 

When moving towards the second singularity from the 
right two additional minima of the potential appear as 
expected for a (supersymmetric) tricritical Ising model. 
Thus for £_|_ 4 - C+,cr we end up with the supersymmetric 
extension of the tricritical Ising model. The only avail¬ 
able scheme is the LPA’ for the cutoff n = 2, for which 
we must tune the anomalous dimension and obtain 


C+(LPA'„ =2 ) = 0.3793822032, 77 = 0.320147. (61) 


The values for all critical parameters are given in Tab. 


XVIII while the critical exponents for even and odd de¬ 


formations are listed in Tables XIX and|XX[ respectively. 



u/(0) 

O'er 

9cr 

96, cr 

7?in 

9out 

V 

LPA’ n= 2 

0.3794 

0.1439 

9.16 

15.36 

0.3201 

0.3201 

1.6653 


TABLE XVIII. Critical parameters for the supersymmetric 
tricritical Ising class in d = 2. 


Shooting 

Sf 

st 

st 

st 

st 

LPA’ b=2 

0.6005 

-0.3129 

-1.5243 

-3.0449 

-4.8593 

SLAC: 

st 

st 

st 

st 

st 

LPA’„ =2 

0.6005 

-0.3129 

-1.5245 

-3.0465 

-4.8710 


TABLE XIX. Critical exponents of the even deformations for 
the on-shell potential of the supersymmetric tricritical Ising 
class in d = 2. 


Shooting: 

st 

S 2 

S 3 

st 

S 5 

LPA’„ =2 

0.8399 

0.1601 

-0.8780 

-2.2468 

-3.9170 

SLAC: 

st 

^2 

^3 

st 

S 3 

lpa+ =2 

0.8399 

0.1601 

-0.8780 

-2.2473 

-3.9213 


TABLE XX. Critical exponents of the odd deformations for 
the on-shell potential of the supersymmetric tricritical Ising 
class in d = 2. 


As for the supersymmetric extension of the Ising class 
in Subsect. |III C[ two critical exponents can already be 
predicted with the help of the results of Appendix [A] In 
the present case the superpotential has the opposite par¬ 
ity as for the Ising class and both deformations discussed 
in the Appendix appear as odd deformations of both the 
potential and the superpotential. The new relations in 
fact determine the first two exponents in Tab. |XX| and 
take the values 

= o*=\, ( 62 ) 


which hold true within our numerical precision. The first 
relation is the equivalent of the superscaling relation (58); 


however, here it does not determine the first critical ex¬ 
ponent of the even deformations and thus does not give 
a direct relation between v and 77 . 


IV. CRITICAL DIMENSIONALITIES 


The two flow equations and ( |41[ ) admit an analytic 
continuation to arbitrary fractional dimensions d > 2, 
leading to the question of whether it is possible to gener¬ 
alize the multicritical models to spaces with fractal prop¬ 
erties. This possibility has been investigated at length 
for the simple scalar model © in H9, whose results we 
partly reproduce below for completeness. However, the 
same is not true for the continuation of the WZ-model 
with one supercharge. One important physical motiva¬ 
tion to continue the system to fractional dimensions is 
that it was observed in [40j that there is a correspondence 
between the scalar models in fractional dimensions with 
local (short-range) interactions, and the same models in 
two or three dimensions with non-local (long-range) in¬ 
teractions. It is thus an intriguing possibility that the 
analytic continuation of our Wess-Zumino model to frac¬ 
tional dimensions plays a role in describing long-range 
effects in the phenomenon of emergent 2 + 1 supersymme¬ 
try in topological superfluid phases [41] . The interesting 
question to answer now is the following: What are the 
critical dimensions at which new critical models appear 
when continuously following d from 4 (where very prob¬ 
ably there is only the Gaussian model) to 2 (where there 
are infinitely many models)? 

The analytic continuations of the WZ-model with more 
than one supersymmetry had already been studied in the 
literature in )42j with the purpose of interpolating be¬ 
tween the a- and E-theorem, and in |43l where the con¬ 
formal bootstrap of the model was considered. The cor¬ 
responding models all contain (on-shell) two real scalar 
fields, have two supersymmetries in 3 dimensions and re¬ 
late to superconformal Af = 2 models in 2 dimensions. 
This family of models is different from those considered in 
the present work. Our models have (on-shell) one scalar 
field and one supersymmetry in 3 dimensions and relate 
to superconformal Af = 1 models in 2 dimensions. 

Motivated by both the Landau-Ginsburg theory and 
the results of the local potential approximation, we as¬ 
sume that the local operators in the effective potential 
are the main actors when it comes to developing a new 
criticality class. A new (upper) critical dimension ap¬ 
pears every time a new local operator becomes marginal. 
To show this, let us first concentrate on the scalar field 
model and equation (15). The set of simple local oper¬ 
ators that satisfy the boundary conditions is exhausted 
by the powers ip 2n for 77, £ N. Each of these local opera¬ 
tors becomes marginal when the corresponding coupling 
is canonically dimensionless; the dimensionality at which 
this happens is thus critical. In the language of (15) the 
most appropriate way to find the critical dimension is by 
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observing that an operator is marginal if 

Sv,r,[VyV'-,<p] | v( ^ )=v , 2n = 0, (63) 

which can be explicitly solved for d to obtain the critical 
dimensions 


d n — 


n 


2 n 


n — 1 n — 1 


V- 


(64) 


The series of critical dimensions accumulates at d = 2 for 
n —> oo when 77 = 0 . 

The plots of the end points <p a of the integration of 
( |15| ) for the first few critical dimensions are depicted in 
Fig- i Each curve can be deformed continuously into 
the others by changing d in (15). It is easy to see that, 


by lowering the dimension, new singularities appear in 
the plot corresponding to new critical models. The way 
in which new singularities appear is as follows: Let us 
start from the upper critical dimension o ?2 = 4 — 277 of 
the ip 4 model where a new singularity appears to the left 
of a = 0. It corresponds to the Wilson-Fisher fixed point 
of the Ising class which we already studied in Subsect. 
|II C[ For d = c ?2 — e with infinitesimal e this singularity is 
arbitrarily close to the Gaussian fixed point solution that 
is generally considered within perturbation theory. If we 
follow d further down to d 3 = 3 — 3/277 the singularity 
has moved left to a finite value of a and just below d 3 
a new singularity appears to the right of rr = 0 corre¬ 
sponding to the tricritical Ising model. Again, this sin¬ 
gularity can be followed for decreasing values of d while 
it moves further to the right. The picture that we just 
described can be iterated further: Each time a new criti¬ 
cal dimension is crossed a new singularity emanates from 
cr = 0. The singularities appear with an alternating sign 
of a = 7/'(0), and this explains the ordering of the mod¬ 
els that we encountered, for example, in Fig. [4j In each 



a 


FIG. 8 . Plot of the function Cp a for the first few critical 
dimensions d n . The plot was obtained using the value 77 = 
1 / 10 . 


interval d n+ 1 < d < d n a new operator ip 2n enters from 
its upper critical dimension defined through (63); thus, it 


is able to satisfy the condition which in turn implies 
a new branching of Fig. [ 8 ] through the condition (18 1 . 

Now we are ready to extend the above discussion to the 
supersymmetric models. Recalling the form of the super¬ 
potential in (351, a new critical dimension is expected to 
occur when a new operator becomes marginal, for 

which the symmetry of the model requires that m £ N. 
Using the notation introduced in (41), we thus seek the 
dimensionality for which 


S WlV [w,w']ip\\ w{ip)=ipm+1 = 0 , (65) 

which can be solved for d to obtain the new critical di¬ 
mensions 


jsusy _ 


2 TO 


m ■ 


m—1 m — 1 


V- 


( 66 ) 


Again, the critical dimensions define a series which ac¬ 
cumulate at d = 2 for m —1 00 when 77 = 0. The two 


sets of critical dimensions (64) and ( 66 ) strictly coincide 


for 77 = 0 , while the difference is related to the fact that 
the on-shell superpotential defined in (37) is subject to 


a multiplicative wavefunction renormalization which has 
no analogue in the scalar field theory. 


Now we can follow the solution of (48) with boundary 


conditions (551 or (541 as a function of the dimension 
when lowering d from d 2 Usy = 4 — 277 to d^ sy = 2 — 77 . 
The result is depicted in Fig.s [9| and [TO] which show the 
function (p for even and odd boundary conditions for the 
superpotential. We see that the dimensions d™ sy for m 



(- 


FIG. 9. Plot of the function ((>q_ for the first few critical 
dimensions d^ sy . The plot was obtained using the value 77 = 
1 / 10 . 

odd are critical for the plot in Fig. [9} while those for m 
even are critical for the plot in Fig. fl0| In both cases, by 
lowering d a new singularity is produced in the plot when 
a new critical dimension is crossed. The alternation of 
the singularities between the Fig.s |9| and [10| is analogous 
to the one that occurs between the left and right sides 
of Fig. [ 8 ] and was expected on the basis of the relation 
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FIG. 10. Plot of the function (p^ + for the first few critical 
dimensions d^ sy . The plot was obtained using the value t/ = 
1 / 10 . 


(56) between the parameters er and (± • Since in our ap¬ 
proximation only the operators of the form Lp n (ipij)) 171 can 
determine the occurrence of new critical dimensions, in 
the limit of zero anomalous dimension our result (66) co¬ 
incides with the set of critical dimensions given in [44j 
where a simple Yukawa model is considered. For this 
reason, it is tempting to speculate that the Yukawa and 
supersymmetric universality classes are the same. This, 
however, is true only to some extent, as the allowed de¬ 
formations of the Wess-Zumino model are constrained by 
supersymmetry if compared to those of a Yukawa system, 
therefore changing the spectrum of the theory (see {45j 
for comments on this model from the point of view of the 
conformal bootstrap). 


V. DISCUSSION 

We studied the functional renormalization group flow 
of simple scalar field theories and their Af = 1 super- 
symmetric extensions in various local potential approx¬ 
imation schemes using a scaling solutions approach, pi¬ 
oneered in [22HB1] . in order to determine the stationary 
solutions of the flow. Our approach to the scalar model 
extends the work in [16] , while the results for the super- 
symmetric Wess-Zumino model are new. Both families of 
models exhibit the expected universality classes in terms 
of stationary solutions of the flow. In three dimensions 
we confirm the presence of both the Ising class and its 
supersymmetric counterpart, while in two dimensions in¬ 
finitely many stationary solutions can be observed which 
we identify with the minimal models of CFT and super- 
CFT. It is rather exceptional that the infinite tower of 
minimal models is visible within such a simple truncation 
of the renormalization group flow that contains only a lo¬ 
cal potential and an anomalous dimension for the field. 

At various stages, we offered a new understanding of 
the mechanism with which new critical solutions appear 


as stationary solutions of the flow of the local potential. 
The key idea is tightly connected to (and best visual¬ 
ized through) the Ginzburg-Landau theory: New critical 
models appear whenever new local operators turn into 
relevant deformations of the spectrum of the theory. We 
explicitly observe that the appearance of each model is 
governed by the canonical scaling of a specific operator 
which can be pinpointed, for example, by analytically 
continuing the renormalization group flow between two 
and three dimensions. The critical models can be identi¬ 
fied as discontinuities in the plot of the domain of exis¬ 
tence of all possible models using both scaling arguments 
and the Landau description that counts the number of 
nodes of the solutions [lj. 

For the sake of the exposition, we concentrated our 
attention on the Ising and tricritical Ising classes, and 
their Af = 1 supersymmetric counterparts. Most impor¬ 
tantly, we introduced two new methods to compute the 
relevant deformations of these models: The shooting and 
the pseudo-spectral SLAC methods. While the shooting 
method has already been applied in [21] to compute only 
the critical exponent of the correlation length, we showed 
that it can be used quite effectively to obtain many low- 
lying exponents. One important ingredient of the present 
work is the application of the SLAC derivative method to 
the computation of the spectrum. This latter approach 
can in fact easily produce a large number of critical expo¬ 
nents with one single calculation, effectively beating all 
other computationally slower algorithms and showing a 
mild dependence on the boundary conditions. The only 
limitation of the SLAC derivative method is in the size 
of the numerically determined critical solution, which is 
finite in our implementation. We thus expect the method 
to work even better when global solutions are available, 
such as in m- 

The functional renormalization group method is an ef¬ 
ficient tool to analytically continue the flow in parameters 
that are fixed in most alternative approaches, for exam¬ 
ple lattice simulations, such as the dimensionality of the 
system. We thus provided an analytic continuation of 
the AT = 1 Wess-Zumino model between two and three 
dimensions that, to the best of our knowledge, is new. 
This continuation can continuously follow the stationary 
solutions of the local potential flow from three to two 
dimensions, showing explicitly what the critical dimen¬ 
sions are at which new models appear. Without wave 
function renormalization the critical dimensions of the 
supersymmetric Wess-Zumino model agree with those of 
the Yukawa model with the same number of on-shell de¬ 
grees of freedom calculated in [33]. It is an interesting 
open question whether this remains true beyond the LPA 
truncation and, in particular, with non-zero anomalous 
dimensions. 

A particular feature of our work is that we can give a 
self-consistent description of any critical model within a 
very simple local potential approximation of the renor¬ 
malization group flow. We provided an extensive and 
detailed numerical analysis of Ising and tricritical Ising 
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models. Our numerical results are in good agreement 
with those obtained by dedicated Monte Carlo simula¬ 
tions and by CFT for the critical exponent of the correla¬ 
tion length but typically fail in determining the anoma¬ 
lous dimension rj accurately. A better determination of 
the anomalous dimension within the functional approach 
is crucial for improving the accuracy of the method. It 
is natural to compare our results for the simple scalar 
field theory with Morris’ contributions to the derivative 
expansion of the renormalization group |22l 124] , in which 
the local potential approximation is complemented by 
a field-dependent wavefunction renormalization and for 
which an additional scaling symmetry (due to the choice 
of a powerlaw cutoff) allows to determine the anomalous 
dimension from the boundary conditions of the flow. We 
find a qualitative agreement with the results of [22] [ 23 ] for 
both the Ising and tricritical Ising classes. As expected, 
the agreement becomes increasingly worse when the field 
dependent wavefunction renormalization of [22] 124] devi¬ 
ates from the value one. The lesson to be learned here is 
that, while the local potential approximation works ex¬ 
traordinarily well in showing that the critical models are 
visible in our functional renormalization group approach, 
better numerical estimates of the critical parameter can 
only be achieved by enlarging the local potential trunca¬ 
tion to include further operators. Therefore, an extension 
of our analysis to truncations of the effective action which 
account for a larger space of operators in functional form, 
such as a field dependent wavefunction renormalization, 
represents the most important and compelling prospect 
of our work. 

Another natural comparison for our work would be 
the recent developments on the conformal bootstrap ap¬ 
proach [Ml HZ3- The latter uses Ward-identities of con¬ 
formal field theory to constrain the operator algebra and 
provide accurate predictions for the critical exponents of 
the Ising model. In our approach, however, conformal 
invariance is an outcome of the requirement of scale in¬ 
variance by spotting stationary solutions of the renormal¬ 
ization group flow and applying the Landau classification 
to the solutions. Thus it is tempting to speculate that in 
a functional renormalization group approach it might not 
be necessary to use the full power of conformal field the¬ 
ory, but simply enhance the correlators appropriately to 
correctly capture the long range correlations in the vicin¬ 
ity of a critical point. Such an extension of the functional 
truncation might roughly follow the lines of [55], or be 
dictated by the desire of maintaining two-loops univer¬ 
sality when continuing the results to four dimensions [38] 
(see also [49] for a conformal field theory perspective). 
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Appendix A: Magnetic and derivative fluctuations 


Let us consider a general renormalization group evolu¬ 
tion equation for a scale dependent “potential” h(p) of 
the form 


kd k h(p) = S h [h,h'-,p]+T h [h"], 
Sh[h , h';ip] = -d h h(p) + d v ipti(ip). 


(Al) 


As done for © we introduced a scaling part Si, [h. h'; p] 
that contains some general (possibly almost canonical) 
dimensions d v and dh for the dimensionful counterparts 
of the field ip and the potential itself respectively. We 
assume that the non-trivial part of the flow Th. \h"] only 
depends on h"(ip) and is otherwise arbitrary. Let us also 
introduce a /c-stationary solution h*(p) of (Al) of given 
parity. The linearized equation for the fluctuations 5h(p) 
defined by 


h(p) —► h*(p) + e8h(<p) — 

. k 0 


-9 


(A2) 


takes the form 


(6-d h )8h(p)+d v p8ti(p)+F h [h*"}5ti'(p) = 0, (A3) 

in which the prime denotes the derivative with respect to 
the argument. 

For any h* (p) with fixed parity there always exist two 
fluctuation modes with fixed parity and critical expo¬ 
nents that are directly related to d v and dh- The first 
one is the magnetic deformation 5h(p) oc ip. It is easy to 
see by direct substitution that this deformation is nor¬ 
malized at the origin by 8 h'( 0) = 1 and solves (A3) for 
the critical exponent 9 = (dh — d^). 

The other fluctuation is 5h(ip) oc h*'(p). To see this 
and to determine the corresponding critical exponent one 
computes 


^(s h [h*,h*'-,p] + F h [h*"}) =0, 


(A4) 


and compares the resulting expression with the fluctu¬ 
ation equation (A3). The comparison reveals that the 


critical exponent of this second deformation is 9 = d v . 
The deformation has opposite parity compared to the 
fixed point solution h*(p). 

In the specific example of the scalar field theory © 
these two deformations have critical exponents 


d + 2 — r] c? — 2 + 77 

--- and --- 

2 2 


(A5) 


The ratio of the second over the first exponent is the 
well known thermodynamical exponent S. The expo¬ 
nents always show up in our numerical computations. 
Since in the scalar case the fixed point solution h* = v* 
is even, both exponents are contained in the spectrum 
corresponding to the odd fluctuations. Given the simple 
relation of these exponents with the anomalous dimen¬ 
sion they can be used, in principle, to test the quality of 
the numerical algorithm used when studying the fluctu¬ 
ation spectrum. 







19 


Appendix B: The LPA’ matching algorithm 


We describe here the matching algorithm that was ap¬ 
plied for obtaining the LPA’ solutions in greater detail. 
Similarly as in the Introduction, we do this for the three 
dimensional Ising class, but the considerations will apply 
straightforwardly to all other examples. 

As discussed in the main text, the stationary solutions 
of equation depend parametrically on the anomalous 
dimension and the required input value when solving this 
equation was denoted by 77i n . For each value of the pa¬ 
rameter r]i n we can compute a corresponding value of cr cr 
by localizing the appropriate singularity in the plot of the 
value (pfj . For reasonably small values of 0 < r} ln < 1 the 
plot of p a is qualitatively similar to Fig. [I] In Fig. 11 we 
plot as a solid line the curve parametrizing the location 
of CT cr as a function of rj- ln . For our definition of LPA’ 
scheme, we need the anomalous dimension computed at 
the minimum po of the solution, which we dubbed 77 ou t) 
to coincide with 77 ^. This is not always the case, as shown 
in Fig. m where the value of 77 ou t is plotted as a dashed 
curve. The mapping cr —>• ? 7 ou t is only defined locally, as 
its definition uses information about the potential at its 
(absolute) minimum, which is not always reached by the 
numerical solution. The intersection of the two curves 
gives, by definition, the condition r]i n = ? 7 ou t, and thus a 
critical solution in the sense of the LPA’. 



FIG. 11. Determination of 77 in the LPA’ of the 3d Ising 
model through the matching condition ? 7 i n = r/ out . The solid 
curve parametrizes the location of cr = a CI as a function of 
77 i n , while the dashed curve is a plot of 77 out of the solution 
corresponding to < 7 cr for the given ? 7 i n . The gray dashed curve 
represents the anomalous dimension of solutions for which 
77 i n = ? 7 out, but it is computed by integrating non critical 
values of a. The dot gives, by construction, a solution which 
is both critical and satisfies the matching condition. 


they make clear that any numerical error in a cr induces 
an error of 77 which is at least one order of magnitude 
bigger. Among the systems studied in this paper, the 
3 d Ising model is the one that shows the smallest error 
propagation from tr cr to 77 . In other cases the amplifica¬ 
tion is by two or three order of magnitudes. Since tr cr 
is the most important parameter as it controls the pre¬ 
cision with which all other parameters are determined, 
we decided to determine it with precision of about 10“ 10 
or higher. The anomalous dimension 77 is then typically 
determined at a precision in the range 10~ 7 — 10 -9 de¬ 
pending on the model considered. 

The Newton-type algorithm for the matching works 
as follows. Once a sufficiently accurate trial value for 
77 ; n ~ 77 (LPA’) is found (for example with precision 10 -3 
as seen in Fig. Ill, this value is used to determine tr cr and 


subsequently 77 ou t- If ?7in and 77 out do not match to the 
desired precision, the value of r]- ln is changed in the direc¬ 
tion of 77 out by a step smaller than the desired precision 
and the procedure is then iterated. The step size should 
be smaller than the desired precision; otherwise, occa¬ 
sionally the algorithm might get stuck in loops, which 
are caused by the fact that the two curves of Fig. EH are 
almost parallel if plotted with the same scale on both 
axes. This very simple algorithm yields precise estimates 
for cr rr . 


Appendix C: Continuity of the critical exponents 


When one tries to calculate the fluctuation spectrum of 
a given fixed point solution one meets the problem that 
the solution of the fixed point equation is only known up 
to a maximal value p a of the field (cf. the corresponding 
discussion in Subsect. II A). This is due to discretization 


and rounding errors involved in the numerical computa¬ 
tion of the critical values of the initial conditions^] Even 
for initial conditions extremely close to criticality one 
finds p < 00 . 

Both the shooting and spectral-SLAC methods have 
been adapted for (various) boundary conditions at the 
maximal value of the field p a for stationary solutions 
of equation parametrized by a to extract the eigen¬ 
function and eigenvalues of the corresponding fluctuation 
operator. This can be done for non-critical solutions and 


the critical solution. In Fig. 12 we plot the exponents 
corresponding to the parity even deformations of the so¬ 
lutions of (151 for d = 3 with boundary conditions (191. 


The plot includes the exponent Of of the constant de¬ 
formation that corresponds to the scaling of the volume 
operator. It was never listed in the tables since for all 
critical models it is just the dimension of the system. 


The plot for the Ising class in Fig. 12 reveals many 
interesting features. All exponents Of seem to be con¬ 
tinuous functions of the parameter a with cusps at the 


The scales of Fig. 11 are particularly instructive, as 


5 Spectral methods do better in finding global solutions 115l |39i . 












20 



FIG. 12. Dependence on a of the critical exponents 9f 
for i = 0 ,..., 6 from top to bottom of the three dimensional 
Ising class. This dependence is computed using the shooting 
algorithm. 


critical values a = <r C r corresponding to the non-trivial 
fixed point and at a = 0 corresponding to the Gaussian 
fixed point. The critical values of cr are marked with ver¬ 
tical dashed lines in the plot. The critical exponents of 
all deformations are negatively divergent at a = — 1 and 
gradually increase with increasing curvature a to their 
value at the Ising class. Beyond a CI they continuously 
interpolate with the spectrum of the free theory at the 
Gaussian point, where they attain their maximal values, 
and then their values decrease again. Each eigenvalue 
shows a distinct behavior near the critical points. Also 
note that the scaling exponent of the volume operator 8 q 
is not equal to 3 below cr cr , but takes this value at cr cr and 
stays there up to a = 0. Beyond this point it decreases 
again albeit very slowly. All higher exponents approach 
cr cr with a slope that differs between the left side and the 
right side, as is expected from the discontinuous change 
in the domain of existence of the potential (see also Fig. 
[2] which shows almost critical solutions from the left and 
the right of the critical point). 


Appendix D: The SLAC derivative method 


We describe here the spectral-SLAC method for the 
computation of the critical exponents that was applied 
in the present work. Again we use the scalar field theory 
of Sect. 0 as an example, but all considerations can be 
extended to the supersymmetric model of Sect. |III| The 
stability analysis around the solution of the second order 
differential equation (15) yields the linear second order 
differential equation (21) that depends parametrically on 


the critical exponents 9. The critical exponents 8 must 
thus, in general, be fixed by the boundary conditions. If 
a global solution to <nn> is present, the physically mean¬ 
ingful condition needed is that the fluctuations are poly- 
nomially bounded in the limit of large field (32], which in 
turn implies a quantized spectrum of critical exponents 


with an upper bound 9 £ {9i}. 

Our numerical solutions, however, do not extend to the 
infinite field (see for example Fig. [3|; therefore, we have 
to impose boundary conditions at the extremum p of the 
domain of validity, rather than at infinity. The influence 
of these auxiliary boundary conditions can be minimized 
by assuming periodicity, in which case the amplitude of 
the fluctuations can vary freely near p rather than be¬ 
ing zero as enforced by the shooting method. With these 
premises, there is a very fast, precise and efficient algo¬ 
rithm to compute the spectrum of the critical exponents 
9i which is known the spectral-SLAC method, because it 
uses the lattice’s SLAC derivative [25] , 

Let the linear ODE describing some fluctuations Su(p) 
around the critical solution v*(p) be of the form 

(-a(p)d 2 +b(p))Su(p) = -95u{p), \p\ < V, (Dl) 


where we used the fact that any linear second order ODE 
can be cast in a form without first order derivative (for 
the scalar model this transformation is given in (22)). 
The functions a(p) and b(p) generally depend on v*Jjp): 
Let us a ssum e that the function a(p) is positive. 

Eqn. (Dl) can be interpreted as the quantum me¬ 


chanical problem of diagonalizing the Hamiltonian "H = 
—a d 2 + b with eigenvalues —9. Now we discretize the 
interval [—p, p\ on which p is defined by N lattice points 
{<£> 1 ,... ,Pn} = A with constant separation 2 p/N. The 
momentum takes its values in the dual lattice 


A - = { W = |(/-I±*) (D2) 

The lattice’s SLAC derivative is defined as 
N . N 

£ fM = -j= , <D3) 

k '=1 VJV 1=1 


where / is a complex-valued function on A and f(p) is its 
Fourier transform. Differently from many standard lat¬ 
tice’s derivatives, the non-local SLAC derivative is given 
by a matrix whose entries are, in general, non-zero for 
any pair of lattice-points pk and pk' ■ 

With respect to periodic boundary conditions the sec¬ 
ond derivative based on the SLAC-derivative reads 



N 2 -l 

3 




(- l) k ~ k ' 


COt jj{k — k') 

sin jf(k — k') 


(D4) 


with an odd number N of lattice points. For anti-periodic 
boundary conditions (suitable for the parity-odd fluctu¬ 
ations) one uses the same matrix but an even number N 
of lattice points. This matrix (D4) is used to discretize 
(Dl) in the form 


^ ^ 7~Lkk' ^ 


(D5) 




















21 


The discretized operator T-Lkk' is an TV x N matrix whose 
eigenvalues can be computed using standard numerical 
methods. Up to a sign these are the exponents The 
method yields very accurate eigenvalues if the corre¬ 
sponding eigenfunctions have their main support within 
the domain \—(p, p\ and if the number of lattice points 
N is sufficiently large to resolve the typical variations of 
the eigenfunction. The biggest advantage of the SLAC 
method is the efficient suppression of discretization errors 
and (for sufficiently large p and N) the very accurate re¬ 
sults for the low lying m N /3 eigenvalues of the spectral 
problem with just one diagonalization. The accuracy can 
be estimated by varying both the number of lattice points 
and the discretization interval [—p,p\. 


Appendix E: The extended supersymmetric solution 


The critical supersymmetric solutions obtained in Sub¬ 
sections |IIIB[ |III C| and |III D| have a limited domain of 
existence, but can be extended to global ones by integrat¬ 
ing through the point p defined by T' w „ = 0 in (51) and 


(|52 ). When approaching the critical point the numerator 
of (481), that represents the scaling of w', fulfills 


^ S w ' ,T] 


0 when p 




(El) 


thus balancing the corresponding zero of the denominator 
and giving rise to a finite w"'(p) in the limit. 

Our strategy for numerically integrating beyond p is 
to perform a linear Taylor expansion of w'(p) and w"(ip) 
for which S w ' >rl and J-' WtV change sign simultaneously at 
the critical point. This can be achieved by choosing the 
almost critical initial condition > C±,cr> for which the 
solution will terminate when T' w „ ~ 0, while S w <is 
still finite. We illustrate in Fig. [13] the behavior of the 


functions T' w and S w p v in the vicinity of (p for the two 
dimensional supersymmetric Ising class found in Subsect. 


Ill C within the LPA’ scheme. The plot shows that the 


two functions almost reach zero at the same point, be¬ 
fore deviating because of the numerical uncertainty in 
the determination of the critical value of £_. The neces¬ 
sary condition for the existence of a finite and non-zero 
w'"(ip ) is that the two functions intersect the <p-axis at 
the same point p. Thus, the crucial point is to accurately 
identify a small intervall in Fig. |13| where T' w v and S v /. rl 
are almost zero but still depend linearly on p. Then 
we perform a linear interpolation that ensures that both 
functions change sign and check that this happens at the 
same point p. 

Once the linear interpolation beyond p is made, the 
numerical integration of (48) yields the solution for all 
values of the field. For large fields the numerator and 
denominator slowly approach zero thus realizing the pre¬ 
dicted scaling limit in (50). In Fig. 14 we show the plot 
of the on-shell effective potential extending beyond p for 
the LPA’ supersymmetric Ising class of Subsect. [me 



0.6150 0.6152 0.6154 0.6156 0.6158 


V 

FIG. 13. The behavior of the functions and S w ' iV near 
the removable singularity. On purpose we chose a value of 
away from criticality, so that the numerical discrepancy 
between the two functions is easily seen. 


The extended solution does not encounter further singu¬ 
larities beyond p even though its asymptotic behavior is 
governed by (50), which is very similar to the asymptotic 
behavior of scalar field theory in (17). 



v(tp) 




FIG. 14. The global on-shell effective poten tial v(p) = 
w'(tp ) 2 /2 for the LPA’ solution of Subsect. IIIC The vertical 


dashed line denotes the position of the singularity i p. 


The natural question to ask at this point is th follow¬ 
ing: Does the extension of the critical solution affect the 
spectrum of its fluctuations? It is possible to show, in 
general, that this is not the case. Let us consider the 
linearized equation for the fluctuations Sw'(p) of w'(tp) 
around the critical solution, which is obtained by deform¬ 
ing 

w'(p) -A w'*(p) + eSw'(p) (E2) 


in the renormalization group flow kdkW r (p) which can be 
easily obtained by applying a p derivative to both sides 
of (41). The linearized ODE is of the form 


/(<P)<9£ + 9 {v)dcp + + 6^ 5w'{p) = 0, 


(E3) 
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in which we introduced the functions 

f(v>) = KA W ") - 

9iv) = F'w^'W'iv) + -—> 


(E4) 


which have to be evaluated at the critical solution w'* (p ). 
As we did in (221 for the simp le scalar model, we trans- 

so that the transformed 


form the fluctuations in 
equation does not contain first derivatives with respect 
to the field. This is done by defining the new fluctuations 
5u{p) such that 


5u(ip) = exp(^ J dyj^yj 6 w\(p), (E5) 

where we can choose po > p or po < p for the lower inte¬ 
gration limit. In both cases, the transformed fluctuations 
solve a Schrodinger-type equation of the form 


Since f(p) = 0 and g{p) = (1 — rj)w" (p) /w"' (p) > 0 the 
transformation (E5) is singular at the point ip. Since in 
addition f(p) has positive slope at p it follows that the 
exponential factor in (E5) vanishes when if approaches ip 
from below and from above. This means that Su obeys 
Diriclrlet boundary conditions at p = ip. It follows that 
the spectral problem (E6) has two disjoint sets of solu¬ 
tions, solutions with p < p and solutions with p > p. 
The only admitted fluctuations satisfy either the bound¬ 
ary conditions (54) or (551, which are possible only for the 
transformed fluctuations for which p 0 < p, thus imply¬ 
ing that p 0 = 0 is the most natural choice. We explored 
numerically the spectrum of the global solution plotted 
in Fig. [14] using the SLAC derivative method disc ussed 
in Appendix [D] applied to the eigenvalue problem (E 61 . 
It can be checked that the spectrum does not receive any 
correction from the extension of the domain of existence 
of the fixed-point solution beyond ip. 


(f(p)dl + P(p)+ 6 ) 6 u(p) = 0. (E6) 
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